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Abstract 


The description of a multifluid code with anomalous transport 
coefficients due to. plasma instabilities self consistently followed 
in space and time is given. As an example we present simulations of 
colliding solar wind streamers. The results compare favorably with 
the observations. 









I. Introduction 

Anomalous transport due to plasma microinstabilities has long 
been known to be of extreme importance in understanding the observations 
of the solar wind plasma . 1 Most theoretical approaches have been 
directed towards the understanding of phenomena that can possibly be 
caused by instabilities or plasma processes 2 (ion heating, anomalous 
wave damping, etc.). Such analytic treatments can provide a rough 
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estimate of the observations and may constitute a useful and necessary 
first step, but they will ultimately break do™ because many different 
plasma effects are acting at once and these effects are not evolving 
in a static spatially uniform plasma as is usually assumed in these 
approaches. The macroscopic properties and plasma microprocesses evolve 
together, feeding back into one. another, and a self-consistent model 
is necessary if the highly non-linear interplay of the interactions is 
to be described properly. 

An increased understanding of the non-linear theory of instabilities 
and recent developments in computational techniques at NRL have made im- 
possible to develop codes that can perform such tasks. These have been 
applied successfully to C.T.R. and laboratory shock problems . 3 In 
this note we describe such a model, which has considerable flexibility 
for application to solar wind problems. As an example we shall simulate 
the interaction of colliding solar wind streams. 

II. - Model ' 

The model consists of a set of multifluid equations, to be 
integrated numerically, which self-cons is tent ly include the effects of 
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wave-particle interactions through anomalous transport terms. These 
terms depend on the instabilities present and evolve in time and space 
as the macroscopic plasma parameters evolve. Although two dimensional . 
models are presently available at NRL, we restrict ourselves for sim- 
plicity to a ID set of equations in Cartesian geometry. Each fluid 
(j) is characterized by its drift velocity u.. its temperature T., 
and its density n. which are functions of space and time. In addition, 
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we follow the evolution of the magnetic field B and the electric field 
E. The coupled set of equations can be written: 
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, equations (1-8) are simply 


a multifluid description of the system coupled collisionally via 
V ja' T jo' given b y the usual Spifczer values. In equation (?) one can 
insert any sources or sinks of energy the particular model may require, 
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as well as thermal conduction. The terms 

represent the effects of anomalous (collective )'proccsses and the 
prescription employed in the code should be derived from independent 
non-linear calculations. We have performed such work for situations 
where electrons and various ion species stream through each other. 

The results have been reported in Ref. 4. In this code we have included 
effects of the magnetized ion-ion two stream instability, the modified 
two-stream instability and the beam cyclotron instability. For these 
instabilities, the momentum transfer terms take the form 
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and the heating terras take the form: 
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where the subscript i refers to the ith ion species, 

j to the j th ion species and e to the electrons. P, p is the total mass 
density, 0^ the lower hybrid frequency, <« is the electron plasma 
frequency, and Fi and F^pj/p ) are factors of order unity. 

These instabilities operate when local conditions are f&vorable. 
The criteria for these instabilities are as follows: 
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Most of the ion heating in the simulations to be described will be. 
due to the ion-ion instability. 
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III. Macro structure of Colliding Streamers 

Using the above model; we proceeded to simulate the regions in 
the interplanetary plasma where fast plasma streams collide with slower 
ones. 5 We employed a three fluid model; one fluid representing the 
fast streaming ions, one the slow ions, and one the electrons. The 
simulation was performed in the reference frame of the slower moving 
plasma, and the initial configuration is shown in Figure 1. Figures 2 
and 5 show the subsequent evolution in time and space, and the agree- 
ment with reported measurements is remarkable. In order to demonstrate 
the effects of the anomalous terms, we equated them to zero in an 
otherwise identical simulation, and the results, shown in Figure 4, 
and are in total disagreement with the observed ion heating. The 
densities, velocities and temperatures plotted in Figures 1-4 are total 
quantities, summed over both ion species. 

In order to demonstrate the mass-dependent differential ion heating 
caused by the instabilities, we performed simulations with proton and 
alpha particle "fast" streams interacting with a proton background. 
Temperatures of the individual species are shown in Figure 5 for the 
proton-proton case, and in Figure 6 for the alpha-proton case. In 
agreement with measurements, T^/T^ ~ 4. 

IV. Micro structure after the Interaction 

The multifluid approach, while it provides a measure of the 
theriualization and of the relative streaming and thermal energies, 
cannot give a description of the velocity distribution function. How- 
ever, since we can find the dominant thermaj-ization mechanism from the 
multifluid calculation, we can proceed to simulate that particular 
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instability in a particle code 6 to find the ion velocity distribution 
function after thermal i z at ion. Figure 7 shows the ion distribution 
before and after thermal izafcion, In the latter case, the distribution 
function is similar to the one measured by Feldman, 7 Cue should note 
that this distribution is no longer unstable, but is marginally stable, 
and will persist over a collisional time scale, the ions streaming 
down the field lines in a one- fluid fashion. What Feldman sees is the 
result of a short wavelength (kR^ « 1) instability which occurred 
at the colliding streamer region and subsequently drifted away. 

Conclusions 

* 

We believe that the multifluid approach outlined here is, at 
this time, the most powerful tool with which a comparison of the 
observations of the solar wind plasma and the various theories can 
be accomplished. The model is flexible enough such that any addi- 
tional processes relevant to the particular problem can be incorpor- 
ated, The example which we presented shows clearly how a combination 
of multifluid and particle codes can explain most of the features 
observed in colliding solar wind streamers, and it demonstrates how 
erroneous conclusions will he drawn if the plasma processes are 
neglected. ^ 
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Figure Captions 


Fig. 1: Magnetic field profiles and total ion temperature, velocity 

and density profiles at t - 0. Background streamer partly 
interpenetrated by "fast" streamer. A sharp initial density 
cutoff is assumed for the "fast" streamer. An initial ion 
temperature of 1 ev, is assumed. 

Fig. 2 : Profiles at t = 10 seconds, for the case with plasma in- 

stabilities operative. 

Fig- 3 : Profiles at t = 20 seconds, for the case with plasma in- 

stabilities operative. 

Fig. 4: Profiles at t 5 20 seconds, for the case with plasma in- 

stabilities turned off. Ion heating results from simple 
adiabatic compression. 

Fig- 5* Temperature profiles for the individual ion streams and for 
the electrons at t - 5 seconds. Here a "fast" streamer 
consisting of protons counterstreams with a proton back- 
ground. 

Fig. 6 : Temperature profiles at. t =.5 seconds, for a "fast" alpha- 

particle streamer counterstreaming with a proton background. 

Fig, 7 ' ion velocity distribution functions before and after 
thermalizat ion via the magnetized ion-ion two stream 
instability. 
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Fig ure 7. Ion velocity d i str i but i on 

f unct i ons . 



